The power of reduced quantum states 
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In a system of n quantum particles, we define a measure of the degree of irreducible n-way 
correlation, by which we mean the correlation that cannot be accounted for by looking at the states 
of n — 1 particles. In the case of almost all pure states of three qubits we show that there is no such 
correlation: almost every pure state of three qubits is completely determined by its two-particle 
reduced density matrices. 



A fundamental question in quantum information the- 
ory is to understand the different types of correlations 
that quantum states can exhibit. A particular issue for 
a quantum state shared among n parties, is the extent 
to which the correlations between these parties is not at- 
tributable to correlations between groups of fewer than n 
parties. In this letter we introduce a way of characteriz- 
ing this irreducible n-party correlation for general states 
of n parties. Our characterization is based on measuring 
the information in the given quantum state of n parties 
that is not already contained in the set of reduced states 
of 7i — 1 parties. 

These considerations lead us to consider the specific 
case of pure states of three qubits. We find the striking 
result that for almost all such states, there is no more in- 
formation in the full quantum state than is already con- 
tained in the three two-party reduced states. Expressed 
differently, the two-party correlations uniquely determine 
the three-party correlations. 

In order to explain our construction, let us first treat 
the case of states of two parties; the local Hilbert spaces 
may have any dimension. Let the (generally mixed) state 
be pab- We ask how much more information there is in 
Pab than is already contained in the two reduced states 
Pa and pB- We address this question by finding another 
state pab which is the most mixed state, i.e. the state 
of maximum entropy, consistent with the reduced states. 
Thus pab contains all the information in pa and pb but 
no more jlj . A simple calculation using Lagrange multi- 
pliers shows that pab has the form 

Pab = exp(A^ <g> l B + 1 A <S> A s ). (1) 

1a and 1b denote the identity operators on the Hilbert 
spaces of particle A and B respectively. A^ and A# come 
from the Lagrange multipliers and are to be determined 
by the condition that the reduced states of pab are re- 
quired to be pa and ps- We can now solve for the La- 
grange multipliers and find that 

Pab = Pa®Pb- (2) 



In the case that pa and ps do not have full rank, this 
calculation is a little delicate since then the Lagrange 
multipliers as they appear in Eq. ([!]) are formally infi- 
nite. In that case we can restrict the Lagrange multipliers 
to the ranges of pa and ps- Then Eq. ([l]) defines Pab 
only on the subspace in which it is nonzero, but Eq. (|^) 
remains valid. 

The difference S(pab) — S{pab), where S is the von 
Neumann entropy, can be interpreted as the amount of 
information in pab that is not contained in pa and ps- 
In fact S(pab) — S(pab) is equal to the quantum mutual 
information S(pa) + S(pb) — S{pab), which measures the 
degree of correlation in pab ■ (Alternatively, we could use 
any measure of the distance between pab and pab to ex- 
press the degree of correlation ||.) We use the word 
"correlation" rather than entanglement since for mixed 
states, pab will have greater entropy than pab if Pab is 
separable but not of product form. For pure states how- 
ever S(pab) = S(pab) if and only if pab is of product 
form, and in this case the difference S(pab) — S(pab) is, 
except for a factor of two, the standard measure of bipar- 
tite entanglement We also note that for a pure state 
with reduced states pa and pb, there are typically many 
states of two parties having the same reduced states. This 
is in contrast to the case for more parties, as we will see 
below. 

We now turn to the more interesting case of quantum 
states of more than two parties; the local Hilbert spaces 
may again have any dimension. For ease of exposition 
we treat the three-party case explicitly; the extension to 
more parties follows straightforwardly. Consider, then, a 
general three-party state pabc- We ask how much more 
information there is in pabc than is already contained 
in the three reduced states pab, Pbc, Pac- 

Before we analyze this situation we point out that there 
are a number of new issues in the three-party case that do 
not arise in the two-party case. Consider a set of states 
Pab, Pbc, Pac which are supposed to be the reduced 
states of some (possibly mixed) state of three parties. 
These three must certainly satisfy some consistency con- 
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ditions: the reduced state pa can arise from both pab 
and paCi and this puts constraints on these two reduced 
bipartite states. However a set of states satisfying this 
condition (and the analogous ones for each of the other 
parties) may still not correspond to a legitimate state of 
three parties. For consider the following set of reduced 
states which are supposed to be the reduced states of 
some state of three qubits: Pab, Pbc, Pac are all sin- 
glets held between the given pairs, e.g. 



PAB = -j=(\0) A \l) B -\l) A \0) B ) 



(3) 



The reduced states of the individual parties are all the 
maximally mixed state of a qubit and so are consistent 
with each other; however it is easy to convince oneself 
that these putative reduced states are not the reduced 
state of any three-party state of three qubits. 

We now return to the main theme of our discussion. 
We are given a general three-party state pabc- We argue 
that a measure of the irreducible three-party correlations 
in the state is the entropy difference between the state 
itself and the three-party state that has no more informa- 
tion in it than in the reduced states. As in the two-party 
case we may use Lagrange multipliers to find the state 
Pabc which contains only the information in the reduced 
states. If pabc has maximal rank, then pabc is of the 
form 

Pabc = exp(A A s <8> lc + Aac ® Is + A BC ® l^)- (4) 

Here Aab, h-BC and Aac come from the Lagrange mul- 
tipliers and are to be determined by the condition that 
the reduced states of pabc be those of pabc- Unlike 
the case of two parties, we have not been able to calcu- 
late these Lagrange multipliers in closed form, in general. 
Nonetheless, the form of pabc is illuminating. For con- 
sider a completely general state of three parties. It can be 
expanded using a basis of operators composed of tensor 
products of operators spanning each individual Hilbcrt 
space. For example, for three qubits, a general mixed 
state may be written as 

Pabc = i(l®l®l + OHOi <g> 1 <g> 1 + fal <g> <Ti ® 1 
+7*1 ® 1 ® (Ji + Ri 3 (Ji <g> Oj ® 1 + S i jO l ®\®Oj 
+Tijl ® <7j <g> c j + Qijk(Ji ® a-j <8> cr fe J , (5) 

since the set of matrices (1, a x , a y , a z ) is a basis for the 
operators on C 2 . It is not the case that the tensor Q de- 
scribes the three-party correlations (for consider a den- 
sity matrix which is of the form pa®Pb®Pc — it has non- 
zero Q). However the discussion above shows that for 
generic density matrices, a state which has all its infor- 
mation contained in its reduced states, has the property 
that its logarithm has no term of the form qijkO~i<g>o~j'g)0~k- 



In a number of places in the above discussion we have 
noted that the case when the states have non-maximal 
rank may need careful treatment. For example one 
clearly cannot take the logarithm of such a state to deter- 
mine whether its information is contained in its reduced 
states. A particularly important class of states of non- 
maximal rank is the set of pure states. As we will now 
see, this set has surprising properties. 

Let us consider the particular case of a system of three 
qubits. All pure states of this system are equivalent un- 
der local unitary transformations to states of the follow- 
ing form Q : 

|97> = a|000> + fo|001> + c|010> + d|100) + e | 111) . (6) 

The labels within each ket refer to qubits A, B and C in 
that order; in what follows we will continue to identify 
qubits only by the ordering of the labels. We now show 
that almost all of these states have no irreducible three- 
party correlation in the sense developed in this paper. 
That is, we show the following: except when the param- 
eters a, 6, c, d, e have certain special values, the state \rj) 
is the only state (pure or mixed) consistent with its two- 
party reduced states. 

Let uj be a three-qubit density matrix whose two- 
particle reduced states are the same as those of \rj). We 
can think of u> as obtained from a pure state \tf>) of a 
larger system, consisting of the three qubits and an envi- 
ronment E: thus uj = Tie | ip) (ip | ■ To get a constraint on 
the form of \ip), consider the state pab of qubits A and 
B as obtained from \rj): 



pab = \4>o)(4>o\ + \4>i)(4>i 



(7) 



where the unnormalized vectors \(f>o) and \4>i) are 

|0o) = a|Q0) +c|01> +d|10); |&> = 6[00> + e |11) . (8) 

We insist that give this same pab when restricted 
to the pair AB. Because pab is confined to the two- 
dimensional space spanned by \4>o) and \4>i), \ip) must 
have the form 



|V;> = \<h)\Eo) + 1^)1^! 



(9) 



where \Eq) and \E{) are vectors in the state-space of the 
composite system consisting of qubit C and the envi- 
ronment E. Computing the density matrix of AB from 
Eq. (||) and comparing it with Eq. (Q), we see that j^o) 
and \Ei) must be orthonormal. It will be helpful to ex- 
pand \Eq) and \Ei) in terms of states of C and states of 
E: 



\Eo) 



|eoo) + |l)|eoi); IE,) = |0)|e 10 ) + |l)|e n ). 



(10) 



Here the environment states |ey) are a priori not nec- 
essarily either normalized or orthogonal. Combining 
Eqs. (|), @ and @, we can write 
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|V) = (a|00)+c|01)+d|10))(|0)|eoo) + |l>|eoi)) 
+ (6|00) + e|ll))(|0)|e 10 ) + |l)|e u )). (11) 

In order to see what further constraints are imposed on 
\ip) by the requirement that the reduced states agree with 
\rj) for the other pairs, let us consider three specific ele- 
ments of the two-party density matrices. 

(11|pbc|H) : As computed from the state \rj), this 
matrix element has the value |< 



As computed from 



Eq. (O), it has the value |c| 2 (e i|e i) 



(en eu . 



(ll|/9Ac|H) : As computed from I77), this matrix ele- 



ment has the value |e| . As computed from Eq. (11), 
it has the value |e| 2 (en|en) + |<i| 2 (eoi|eoi). Hence for 
generic values of c, d and e, |eoi) = and (en|en) = 1, 
from which it follows that |exo) = and (eoo|eoo) = E 

(01|/?bc|10): As computed from \rj), this matrix el- 
ement has the value be*. As computed from Eq. (fill ) 
(with |eoi) = |eio) = 0), it has the value &c*(eoo|en). 
We conclude, again for generic values of the parameters, 
that |e 00 ) = \eu). 

Inserting these inferences into Eq. (|ll]), we find that 

\ip) = (a|000) + 6|001> + c|010) + d|100) + ejlll»je 00 >- 

(12) 

When we trace out the environment to get the state w, we 
see that we must have u> — \r/)(r/\. That is, the only state 
(pure or mixed) consistent with the two-particle reduced 
states of I rf) is \rj) itself. 

The above treatment deals simply with the generic 
pure state of three qubits. We have found it necessary to 
use a slightly more involved analysis, to be found in the 
Appendix, to identify those special states for which the 
two-party reduced states do not uniquely determine the 
full three-party state. The results in the Appendix show 
that the only states that do not have this generic prop- 
erty are those which are equivalent under local rotations 
to states of the form 



a|000) +6|111). 



(13) 



We can identify these exceptional states as the ones that 
admit a three-particle Schmidt decomposition as studied 
by Peres || . 

The results of this letter clearly raise many questions. 
For example whether the properties that we have found 
for generic pure states of three qubits extend to systems 
of more parties and in higher dimensional Hilbert spaces; 
we intend to return to this in a future publication. Also it 
is interesting to find non-trivial classes of n-party states 
that are determined by their reduced states of fewer than 
n — 1 parties, and to characterize their entanglement 
properties. An example is the family of states 

o|0001) + 6|0010) + c|0100) + d|1000). (14) 

These states are uniquely determined by their two-party 
reduced states. 



Finally, we note that many of the ideas we have put 
forward here also shed light on classical probability dis- 
tributions. For example the idea of characterizing the 
n-party correlations using the information in the (n — 1)- 
party marginal distributions. In the light of our results 
on pure states of three qubits it is intriguing to consider 
the case of probability distributions P(X, Y, Z) of three 
random variables each of which has two values; such a dis- 
tribution arises from local von Neumann measurements 
on states of three qubits. In this case it is not difficult 
to see that generic distributions are by no means deter- 
mined by their marginal distributions. For consider a 
given set of probabilities pijk where pooo is the probabil- 
ity that X = 0, Y = 0,Z = etc. The set of probabilities 
Qijk = Pijk + bi—^y^ 1 ^ has the same two-party marginal 
distributions, where 8 is a constant and e(ijk) is the par- 
ity of the bit string ijk. 
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Union under the project EQUIP (contract IST-1999- 
11063). 

Appendix 

Consider an arbitrary pure state = Ylijk a ijk\ijk) 
of three qubits A, B and C. We give an alternative 
derivation that for generic a,ijk, \rf) is uniquely deter- 
mined by its two-party reduced states and find those 
states for which this is not true. 

We can quickly dispose of the case in which \rj) is the 
product of a single-qubit state and a two-qubit state. 
In that case the two-party reduced states determine 
both factors in the product and therefore determine \rf) 
uniquely. In what follows, we will assume that \rj) does 
not have this product form. 

A general state that agrees with \rj) in its reduced 
states can always be obtained from a pure state \ip) of 
the three qubits plus an environment E. Let us first ask 
what form \ip) must take in order to be consistent with 
the (generally mixed) state of the pair AB derived from 
I rf). By an argument essentially identical to the one lead- 
ing to Eq. (|ll|), we find that must be of the form 



1^ 

ijkl 



aijl\ijk)\ei k ) 



(15) 



Here I takes the values and 1, and the states |e/fc), 
which are states of E alone, satisfy the orthonormality 
condition 

/~2(eik\ei>k) = 5u>- (16) 



Similarly, by considering AC and BC we see that 

\ip) = /^2aiik\ijk)\fij) =Y^ a w\ i i k )\9ii), (17) 

ijkl ijkl 
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with J2j(fij\fi'j) = S w and J2i(9n\gi'i) = Here we 
regard the coefficients a^ as fixed — that is, the state 
\rf) is fixed — and we are looking for environment vectors 
\ e ik), \fij) and \gu) that satisfy the various linear equa- 
tions arising from the fact that the three expressions for 
\ip ) in Eqs. (15) and ( |l7| ) must all be equal. 

It is instructive to write down explicitly, as an exam- 
ple, the two equations arising from (|l|) and (0) that 
involve only the vectors |e 00 ), |e 10 ), |/oo) and |/i ): 

aooo|eoo) + aooi|eio) = aooo|/oo) + aoio|/io) (18) 
aioo|eoo) + aioi|eio) = aioo|/oo) + aiio|/io) • (19) 

Notice that these two equations are linearly independent: 
if they were not, then the state I77) would be factorable 
into a single-qubit state and a two-qubit state, contrary 
to our current assumptions. 

These equations and analogous ones relating \eik) to 
\gu) and \fij) to \gu) can be solved fully, and one finds 
that the general solution for the vectors \eiu) is 

|eoi) = (aonaioi - a m aooi)|u) 

|eio) = (aoooano - ctiooaoio)|v) 
|eoo) = (aioodon + aioiaoio)^) + \w) (20) 
|en) = (aoooam + aooiaiio)|v) + \w) , 

the vectors \v) and \w) being arbitrary. The correspond- 
ing expressions for the / and g vectors can be obtained 
from Eq. ( po| ) by permuting indices; for example, the ex- 
pression for each / vector is obtained from the expression 
for the corresponding e vector by permuting the last two 
indices of every ayfc (without changing the vectors \v) 
and \w)). Thus, once the two vectors \v) and \w) have 
been chosen, the solution is determined. The form of 
the solution shows that at most two dimensions of the 
environment can ever be used. 

Wc have not yet taken into account the orthonormality 
conditions for the environment states. Let us now con- 
sider just Eq. ( [l6| ) which constrains the e vectors. It is 
helpful to rewrite Eq. (|2^) in terms of a new arbitrary 
vector \z) that replaces \w): 



|e i) = a\v) \e w ) = P\v) 
|e 00 ) = \z) |en) = \z) +-y\v) 



(21) 



where a = Oonfflioi - omaooij P = aoooano - aiooaoio 
and 7 = aoooam + ^001^110 — aiooaon — aioiaoio- In 
terms of these parameters, the orthonormality condition 
of Eq. (|l^) is expressed by the following three equations 

(z\z) + \a\ 2 (v\v) = 1 
(z\z) + (\p\ 2 + h\ 2 )(v\v) + l(z\v) + j(v\z) = 1 (22) 
aj(v\v) + P(z\v) + a(v\z) = 0. 

Taking the difference between the first two of these equa- 
tions, and treating separately the real and imaginary 
parts of the third, we obtain three homogeneous linear 



equations for the three real variables (v\v), Re (z\v) and 
lm(z\v). For generic values of a, (3 and 7, these three 
equations are linearly independent, so that the only so- 
lution is \v) = 0. This in turn implies, by Eq. 
that |eoi) = |eio) = and |eoo) = |en) . Thus in this 
generic case only a single dimension of the environment 
is used — that is, the environment is in a pure state — and 
the qubits ABC must be in the given state \rj). 

This conclusion can be avoided only if the determinant 
D of the 3x3 matrix associated with the three homo- 
geneous linear equations vanishes, and the corresponding 
determinants computed from the two other orthonormal- 
ity conditions (for the vectors / and g) are also zero. 
Computing D explicitly, we find that D = if and only 
if (i) \a\ = \P\ and (ii) 7 2 a/3 is real and non-negative. 

Suppose now that I77) is not determined by its two- 
party reduced states, so that the above conditions (i) 
and (ii) must be satisfied. These conditions imply that 
there exists a local rotation on qubit C that will bring 
both a and /? to zero, thus bringing \rj) to the form 



|PA)bB>|0) + |gA>|9fl)|l) 



(23) 



Here \pa) and are (unnormalized) vectors in the 
space of qubit A, and \pg) and \qs) belong to qubit B. 
We now use in a similar way the conditions analogous to 
(i) and (ii) but derived from the orthonormality relations 
for the / vectors. These imply that we can apply to the 
form ( p3| ) a local rotation on qubit B to bring it to the 
form |p^)|0)|0) + 1*7^4) 1 1) 1 1) . Finally, from the conditions 
derived from the g vectors, it follows that we can rotate 
qubit A and arrive at the form a|000) + 6|111). 

We conclude, then, that the only pure three-qubit 
states that might not be uniquely determined by their 
two-particle reduced states are those that are equivalent 
under local rotations to the form given in Eq. ([l3]) . In fact 
it is easy to see that for any state of this form with b^O 
and b ^ 0, there do exist other three-qubit states — e.g. 
mixtures of 1 000) and | 111) — having the same two-party 
reduced states. 
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